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ABSTRACT 

We present a study of the line-of-sight magnetic fields in five large-diameter 
Galactic H II regions. Using the Faraday rotation of background polarized radio 
sources, as well as dust-corrected Ha surface brightness as a probe of electron 
density, we estimated the strength and orientation of the magnetic field along 93 
individual sight-lines through the Hll regions. Each of the Hll regions displayed 
a coherent magnetic field. The magnetic field strength (line-of-sight component) 
in the regions ranges from 2 to 6 fiG, which is similar to the typical magnetic 
field strength in the diffuse interstellar medium. We investigated the relationship 
between magnetic field strength and electron density in the 5 Hll regions. The 
slope of magnetic field vs. density in the low-density regime (0.8 < n e < 30 cm~ 3 ) 
is very slightly above zero. We also calculated the ratio of thermal to magnetic 
pressure, (3 t h, for each data point, which fell in the range 1.01 < (3 t h < 25. Finally, 
we studied the orientation of the magnetic field in the solar neighborhood (d < 
1.1 kpc) using our data from 5 Hll regions along with existing measurements of 
the line-of-sight magnetic field strength from polarized pulsars whose distances 
have been determined from their annual parallax. We identify a net direction for 
the magnetic field in the solar neighborhood, but find no evidence for a preferred 
vertical direction of the magnetic field above or below the Galactic plane. 

Subject headings: ISM: magnetic fields — polarization — radio continuum: ISM - 
H II regions 
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Introduction 



Understanding the interactions between neutral and ionized matter and magnetic 
fields is crucial in testing models of star-formation in the Galaxy. The expectation from 
ideal magnetohydrodynamic flux freezing is that the total magnetic field strength, B, is 
proportional to n K , where n is the gas number density and < k < 1, depending on the 



geometry of bulk gas motions relative to the magnetic field lines (Crutcher 1999). To date 



observational estimates of magnetic field strengths in a large number of Galactic molecular 
clouds suggest that k ~ 0.4 — 0.6 in regions with volume densities between 10 3 — 10 5 cm -3 



(Troland & Heiles 1986; Vallee 1995 Crutcher 2007; Crutcher et al. 2010). This range of k 



probably represents an upper limit, as processes such as fast ambipolar diffusion (ion-neutral 



drift) may also play an important role within dense, turbulent regions (Zweibel 2002). In 
contrast, regions with lower gas densities (in the range of 10 _1 — 10 3 cm" 3 ) appear to have 



magnetic fields of a fairly constant ~5 — 10 /iG (Crutcher 2007). A simple interpretation 



is that, in the absence of strong turbulence, low-density gas in the ISM is free to move 



along magnetic field lines without inciting a change in magnetic field strength (Fiedler & 



Mouschovias 1993 Padoan & Nordlund 1999 Kim & Ostriker 2006). The constancy in 



the strength of the magnetic field at low densities has been predicted by models (Troland 



& Heiles 1986; Fiedler & Mouschovias 1993; Heitsch et al. 2009), but the mechanism by 



which material accumulates without increasing the magnetic flux is not clearly understood. 



Numerical simulations by Heitsch et al. (2007, 2009) have shown that gas flowing parallel to 



magnetic field lines may generate strong turbulent instabilities which, through the depletion 
of magnetic flux by reconnection or magnetic flux dissipation, may actually enhance the 
correlation between B and n rather than destroying it. A major obstacle to understanding 
the process of large-scale star formation is the lack of reliable measurements of magnetic 



fields at relatively low gas densities (Troland & Heiles 1986; Crutcher 2007). Particularly 



valuable are in situ measurements of B and n across a range of gas densities. 
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Measurements of magnetic fields at a range of low gas densities also enable us 
to evaluate the magnetic and thermal pressures, which are useful inputs to numerical 
simulations of galaxy structure. These models often assume that the weight of the gas 



layer in the Galactic disk is counteracted by the total pressure of the disk gas (Piontek & 



Ostriker 2007). Ostriker et al. (2010) presented a theory of star formation in a multi-phase 
interstellar medium (ISM) that is self-regulated by stellar heating. The ratio of total to 
thermal pressure in low-density gas is a vital input to this model, but is poorly constrained 



by observations, and estimates quoted by Ostriker et al. (2010) vary by an order of 



magnitude. We can assist these models by providing measurements of the thermal and 
magnetic pressures in the ISM over a range of densities. 

Calculating the magnetic pressure requires a robust in-situ measurement of magnetic 
field strength, which is a difficult task. An excellent probe of magnetic fields is Faraday 
rotation of extragalactic radio sources, provided that the Faraday rotation due to the 
magnetized plasma in the region under scrutiny can be separated from the Faraday rotation 
due to the remainder of the Galactic interstellar medium along the line of sight. By focusing 
on discrete regions of relatively high electron density, i.e. Hll regions, it is possible to detect 
the Faraday rotation 'signal' due to these localized regions above the 'noise' provided by 
Galactic foreground and background. With this in mind, we here present a study of a small 
sample of large-diameter H II regions, which we use to derive the magnetic field strength 
and (5 t h across a range of gas densities. 

In this experiment, we have used Faraday Rotation of background polarized radio 
sources and previously published maps of Ha surface brightness to measure the magnetic 
field and gas density in 5 Hll regions that lie within 1.1 kpc of the Sun. By focusing on 
discrete regions of relatively high electron density and calibrating against the mean Faraday 
rotation along similar sightlines, it becomes possible to isolate the Faraday rotation due 
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to the H II region. The technique we use is based upon the foundations laid by |Heiles 



Chu (1980) and Heiles et al. (1981), who estimated magnetic field strengths in Hll regions 
using the Faraday rotation of just one or two background polarized sources. With a much 
larger number of background polarized sources (93 in total) to choose from, we are able to 
significantly advance this technique to accurately measure the line-of-sight component of 
the magnetic field strength in 5 discrete ionized volumes of the Galaxy. 

This paper is organised as follows. In Section 2 we introduce our method for measuring 
B\\, n and j3 t h in Hll regions. Section 3 contains a description of our source selection process 
and our model for the Hll regions. Section 4 describes our results, including maps of the 
H II regions, plots of rotation measure vs. Ha surface brightness and plots of B vs. n. In 
Section 5, we compare our work with previous results, examine the relationship between 
magnetic field strength and gas density in our H II regions, and estimate the corresponding 
ratio of magnetic to thermal pressure. We comment on the orientation of the large-scale 
magnetic field in the local Galaxy, both above and below the plane. Our conclusions are 
presented in Section 6. 



2. Experimental Details 

2.1. Probes of magnetic field and electron density in Hll regions 

H II regions are Faraday screens, i.e. the magnetized ionized plasma causes the electric 
field vector of a linearly polarized wave traversing the region to rotate by an angle that 
is proportional to the square of the observed wavelength, A 2 . Rotation Measure, RM, is 
the scaling factor between A 2 and the change in polarization position angle, Xi due to 
Faraday rotation, i.e. Ax = RM A 2 . Quantitatively, the RM of an Hll region is described 
by Equation Q: 
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RM = 0.81 J n e B\\ dl rad m~ 2 (1) 

where n e is the density in cm -3 of free electrons in the Faraday rotating medium, B» is the 
parallel component of the magnetic field in jiG and the integral along a line element dl runs 
from the source to the observer in parsecs. Positive and negative RMs indicate average 
magnetic fields pointed towards and away from the observer, respectively 

The rotation measure of polarized radiation from a background source can be used as 
a probe of the electron-weighted parallel component of the magnetic field strength along a 
single sightline through an Hll region. The most straightforward way to infer B\\ from RM 
is to use a background pulsar. The dispersion measure, DM, of a pulsed electromagnetic 
signal from a pulsar is given by DM = n e dl cm -3 pc, where the integral limits are 
the positions of the observer and the pulsar. The average magnetic field strength along 
the line of sight to the pulsar is simply given by B\\ = ■ However, the measured 
DM represents the sum of free electrons inside the H II region and in other gas along the 
same sightline. Therefore the magnetic field within the Hll region may only be inferred if 
a foreground subtraction can be made for DM using other pulsars at a similar distances 
in the same region of sky that do not intersect the H II region. Similarly, a correction for 
Faraday rotation in the Galaxy is made using RMs of polarized extragalactic sources that 
lie close to, but outside the H II region boundary. Unfortunately, the number of pulsars for 
which distances have been accurately determined is sufficiently low that this situation does 
not commonly arise. 

There is another source of information on rotation measure aside from pulsars — 
namely the increasingly large sample of unresolved polarized extragalactic sources. All 
extragalactic sources lie at essentially the same distance (i.e. at infinity) relative to the 
magneto-ionized material in the Galaxy, so reliable corrections for the Galactic RM can be 
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made using an ensemble of extragalactic sources in a surrounding region of the sky. We use 
the RMs of polarized background extragalactic point sources to probe the electron-weighted 
magnetic field strength of several different paths through H II regions. Our source of data 
on Faraday rotation is a recently published catalog of the RMs of 37,543 polarized radio 
sources ( |Taylor et al.|2009] ) derived from the NRAO VLA Sky Survey (NVSS; |Condon et aL 



1998) at 1.4 GHz. 



2.2. Source Selection 



To identify a suitable source sample for this study, we searched maps of the Ho; sky 
for prominent Hll regions. To do this, we used maps of Ha surface brightness I// Q , with 



6' resolution (Finkbeiner 2003) compiled from the WHAM (Haffner et al. 2003), SHASSA 



(Gaustad et al. 2001) and VTSS (Dennison et al. 1998) surveys. Our search was restricted 



to regions at least 5° away from the Galactic plane, as the RM signal close to the plane is 
dominated by the dense magnetized plasma in the Galaxy and dust extinction maps used to 



correct emission measures are unreliable (Schlegel et al. 1998). We identified 5 Hll regions 



which lie at Galactic latitudes |6| >5°, which have a peak ln a at least 3 times the root 
mean square surface brightness of the surrounding region and have at least 4 background 



RM measurements in the catalogue of Taylor et al. (2009). These were Sh 2-27, Sh 2-264 



Sivan 3, Sh 2-171 (NGC 7822) and Sh 2-220 (The California Nebula). Table [T] shows the 
basic parameters for each region adopted in our analysis, including the number of RMs 



from Taylor et al. (2009) lying behind each region, the diameter of the Hll region, and 



properties of the ionizing star: distance, reddening and spectral type. For each H II region, 
we defined the boundary as the line within the regions of sky shown in Figure [T] for which 
Iffa = 3 (TrmSi where ctrms is the root mean square of I# a outside the boundary (see details 
in Table [21. 
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2.3. Calculating B\\ for our Hll region sample 



In Section |2.1| we stated that RM probes the electron-weighted magnetic field strength 
along the line-of-sight. Therefore, by estimating the electron density we can disentangle 
B\\ from RM in Hll regions. We do this by measuring the Emission Measure, EM which is 
defined as EM = J °° v? e dl cm -6 . There are several observational probes of EM, the most 
common being the n = 3 — 2 Balmer (Ha) recombination line of atomic hydrogen at 6563A. 



In order to determine the EM corresponding to each RM from the Taylor et al. (2009) 



catalogue, we used the Finkbeiner (2003) compilations of I# Q to estimate the electron 



density corresponding to each measured value of RM. 



The EM along a path is related to I# a (Reynolds 1988) in the following way: 



EM = 2.75 



10 4 K 



0.9 



<-Ha 

"r" 



t -6 

e cm pc. 



(2) 



Here, T e is the electron temperature and In a is in rayleighs, where 1 R=10 6 /47r 
photons s _1 cm~ 2 sr _1 . A correction term for dust extinction, e T , where the optical depth 



r = 2.44 x Eb-v (Finkbeiner 2003), is included to extract the intrinsic Ijj a from the 



measured value. This assumes that all the dust is in front of the source, rather than mixed 
in the source or behind the source. Several features in the dust maps towards our five Hll 
regions were also visible in the Ha maps, indicating that much of the dust is in front of 
the Hll regions. We therefore applied the maximum possible correction to I# a , assuming 
that all the dust is situated between the H II region and the Earth. Our dust corrected EMs 
should therefore be regarded as upper limits. 

We adopt the following formalism to allow for possible variations in electron density 
along any line of sight through an H II region. Our model assumes that along a single line of 
sight with depth L pc, the ionized material is composed of discrete clumps. These clumps 



-9- 



occupy a total length /L, where / is the volume filling factor. Where ionized material 
is present (i.e. within a clump), the electron density is n e = n . Outside the clumps we 
assume that n e = 0. 

The electron density in a clump is given by 



/EM o 

n ° = V ~Jl cm ' ^ ' 

and the average density along a sightline is (n) = fn . 
In this case the strength of the parallel component of the magnetic field in the H II region is 



fill = ™ ^_ fxG, (4) 



assuming that n and Bn are uncorrelated along a single sight-line. If fluctuations in the 
electron density and magnetic field are positively (negatively) correlated, then the use of 
Equation Q will lead to an overestimate (underestimate) of the magnetic field strength 
(see 



Beck et al. 2003). We will discuss this point further in Section 5.3. 



To evaluate Equation Q, RM and EM are determined from observations. The total 
path through the Hll region, L, can be estimated from the geometry of the region and / 
may be assumed to within some reasonable physical limits. Therefore, using Equations (|3| 
and (El) we can estimate n and Bu at the position of each RM measurement. 



Once these quantities are known we can calculate the magnetic pressure P mag , 
thermal pressure Pther an d hence the plasma beta ((3 t h= Pther/Pmag) for each data point as 
follows. Using the statistical approximation B to t ~ \/3 Bn and the observational estimate 



(B uni f orm / B tot ) < \/2 (Beck 2001) we first estimate the total magnetic field strength 



Btot ~ V^Bn from the measured magnetic field strength, Bn for each data point. The 
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magnetic pressure P m ag= Bf ot /8TT and the thermal pressure, Pther = 2n kT e can then 
be calculated, where k is Boltzmann's constant and assuming T e =7000 K and a 100% 
ionization fraction. Finally (3th can be determined for each data point. 



2.4. Background pulsar search 



In Section 2.1 we explained how measurements of the RM and DM of pulsars can 



directly yield B\\ along a line of sight. In order to identify possible pulsars intersecting our 
sample of H II regions we carried out a search of the ATNF pulsar database^] ( [Manchester 



et al. 2005). Of our sample of five Hll regions, Sh 2-27 has three pulsars and Sh 2-264 



has one pulsar that overlaps the H II region. DM-derived distances from the pulsar catalog 



are unreliable, as they rely on models of electron density in the Galaxy (Taylor & Cordes 



1993 Cordes & Lazio 2002) that do not include the Hll regions in our sample. In order to 



ensure accurate distance estimates, thus selecting only pulsars that lie behind H II regions, 
we selected only the pulsars with distances derived from measurements of annual parallax. 



This left a single pulsar, PSR J1643-1224 (Lorimer et al. 1995), which lies behind Sh 2-27. 



We perform a consistency check on our results for Sh 2-27 in Section |3j using the measured 
properties of this pulsar. 



3. Results 



Figure [T] shows the observed Ha surface brightness (grayscale) and rotation measures 
(circles) for each of our Hll regions, displayed in a stereographic projection. Blue and red 



1 On-line catalogue at 
November 2009 



http : //www. atnf . csiro . au/research/pulsar/psrcat 



v.1.38, 
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circles indicate negative and positive RMs, respectively. The position of PSR J1643-1224 
behind Sh 2-27 is denoted by a filled green circle. Ha maps were smoothed to a resolution of 
1° (the resolution of the WHAM beam). For each background polarized source we measured 
the Ina of the Ha pixel closest to the NVSS source. In the case of two overlapping pixels 
the source was assigned the mean l Ha of the two. 

It appears that the magnitude of the RM of polarized sources is greater inside the H II 
region boundary than outside the boundary, particularly for Sh 2-27 and Sh 2-264 with the 
highest Ha surface brightnesses. In order to inspect this relationship more clearly, we made 
plots of RM vs. lHa for each region. Figure [2] shows these plots of the observed RM vs. 
ln a for polarized background sources in each region of sky shown in Figure [1} Black circles 
indicate NVSS sources that lie outside the boundary of the Hll region. Red and blue circles 



denote sources from Taylor et al. (2009) with positive and negative RMs, respectively, that 



lie behind the Hll region. 

The plots in Figure [2] show unambiguously that, for polarized sources whose radiation 
has passed through an Hll region, there is a positive correlation between I# a and |RM|. 
This is most clearly evident for Sh 2-27 and Sh 2-264, which have a large number of 
background polarized radio sources. The fact that RM has a consistent sign (positive or 
negative) within each region indicates that the magnetic field is largely coherent. Such a 
large-scale coherent magnetic field in the ionised gas of H II regions has been predicted by 



the MHD simulations of Arthur et al. (2011). In addition to this coherent signal there is a 



scatter in RM ~ ± 100 rad m 2 , which is likely related to the random component of the 
magnetic field. We defer our discussion of random magnetic fields to a future paper. 

For Sh 2-220 (the California Nebula), inspection of Figure [IJe) initially suggests that 
the Hll region is introducing a negative RM to the extragalactic sources. However, in 
Figure |2^e) we also see four positive RMs between 60—90 rad m -2 that apparently lie 
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within the boundary of the H II region. Closer investigation reveals that these points lie on 
the northern edge of the Sh 2-220 H II region, which has a very sharp drop-off in Ha surface 



brightness due to the interaction between the Hll region and a molecular cloud (Elmegreen 
fc Elmegreen||1978 ). The diameter of the WHAM beam is 1° and any polarized source along 
the edge of the H II region is assigned a mean value of for the whole WHAM beam. We 
therefore assume that the Ha surface brightness plotted on Figure [T] for these four sources 
is artificially high due to beam averaging effects. Future, high-resolution Ha observations 
could confirm this interpretation. 

Next, we present the calculations of electron density, no and the magnetic field parallel 
to the line-of-sight, Bu. The measured value of ln a includes emission from the Hll region 
plus contributions from the foreground and background warm ionized material in the 
Galaxy. In order to account for this, we calculated the mean Ha surface brightness from 
a representative region of sky surrounding the H II region boundary and subtracted this 
value from all Ha surface brightness measurements within the boundary of the H II region. 
Similarly, the extrinsic RM (or in some cases, a systematic RM gradient across the field) 
was subtracted from all RM values within the Hll region boundary. Details of the ln a 
boundary conditions and magnitudes of the extrinsic RM and corrections in the Hll 
regions are given in Table [2] 

After correcting the measured values of RM and ln a for foreground and background 
contributions, we converted I# Q to EM using Equation (j2J), assuming that T e = 7000 K 



(Reynolds 1988 Madsen et al. 2006) and adopting values of E#-y from the literature 
(listed in Table [T]). We then calculated the electron density, no, of the Hll region along the 
sightline to each background polarized source. We defined a central position for each Hll 
region and calculated the path- length through the region for each background source, using 
the equation L = 2[R 2 — D 2 sin 2 {6)] 1 / 2 where R is the radius of the sphere, D is the distance 
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from the observer to the center of the sphere and 9 is the angle between the line that goes 
from the observer to the center and the line that goes from the observer to the background 
source in question. The electron density no derived from each Ha measurements was 
then calculated using Equation (j3|. We adopted a volume filling factor / = 0.1, which is 
representative of the filling factors estimated for a number of H II regions of this size in the 



literature (Kassim et al. 1989 Herter et al. 1982). Finally, we calculated Bn for the Hll 



region towards each polarized source using Equation Q. 



Figure [3] shows a plot of log n against log B\\ for our 5 Hll regions. The black line 
represents equilibrium between magnetic and thermal pressures, which is discussed in 



Section 4.3. There is a significant scatter of both no and Bu within each region. This scatter 
in Bu is possibly related to turbulence, as the magnitude of Bn changes with the direction of 
the magnetic field vector with respect to the line of sight. Some of this scatter may also be 
due to the uncertainty in Eb-v along individual lines of sight, which affects the calculation 
of Bu. Another contributing factor to the scatter in the data may be the uncertainty in 
the electron filling factor for each sightline. In Figure [3] (bottom right) we have indicated 
with an arrow the direction and distance by which each data point will move if the filling 
factor changes from / = 0.1 (plotted) to values of / = 0.04 and / = 0.25. This range in 
/ is representative of filling factors measured in Hll regions. The standard deviation of 
the optical depth, a T , for each Hll region is displayed as an arrow on the left-hand side of 
Figure |3j Colored arrows (left) indicate the distance and direction each data point would 



move were the optical depth reduced by one standard deviation. 



Table |3j lists the derived values of the diameter, B\\ and no of each Hll region, where 
bars indicate median values for each Hll region. The estimated median electron densities 
for each region (approximately equal to the density of hydrogen atoms, n(H), assuming a 
fully ionized medium) have values 1.5 < no < 14.1 cm -3 . The line-of-sight magnetic field 
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strengths in our Hll regions range from —6.3 to +2.2 /zG, similar in magnitude to estimates 



of the magnetic field strength in the diffuse ISM using other methods (Heiles & Crutcher 



2005 Crutcher 2007). This suggests that the magnetic field in these large-diameter Hll 



regions may be the ambient magnetic field of the Galaxy that is 'lit up' by virtue of the 
high electron density of the Hll region. Estimates of the uncertainties in B\\ and n are 
quoted, although the derived magnetic field strengths may be over- or underestimated 
by approximately a factor of three, due to possible correlation of fluctuations in electron 



density and magnetic fields in turbulent media (Beck et al. 2003) 



Following our analysis, we carried out a consistency check on our results by using the 
RM and DM of PSR J1643-1224 to independently measure the filling factor and magnetic 



field strength of Sh 2-27. Pulsar J1643-1224 was discovered by Lorimer et al. (1995) and 



lies at a distance (derived by annual parallax) of 455 ± 144 pc (Verbiest et al. 2009). The 
RM of PSR J1643-1224 is -263± 15 rad m" 2 ((Manchester & Han||2004l) and its DM is 62.4 



cm pc (Toscano et al. 1999). We removed the foreground contributions to the RM and 



DM due to the magnetoionic medium of the Galaxy in the direction of PSR J1643-1224 by 
subtracting the RM = +6.5± 10 rad m~ 2 (|Manchester||1974|) and DM = 10.7 cm" 6 (iHobbs 



et al. 2004) of PSR B1604-00 from the measured values for PSR J1643-1224. We used PSR 



B1604-00 for this correction as it lies at a similar distance to PSR J1643-1224 (590 ± 148 



pc), derived by the electron density model of Taylor h Cordes (1993) and has co-ordinates 



close to, but outside the boundary of Sh 2-27. 



Putting the foreground-corrected RM and DM into the relationship B\\ 
find that B» = —6.4 /iG. This is consistent with the value B» 



RM 



we 



0.81DM' 

-6.1 ± 2.8 fiG that we 

derived from polarized extragalactic sources, within the uncertainties. Ha composite maps 



published by Finkbeiner (2003) allowed us to calculate a residual ln a (after subtracting the 
mean background of 0.9 R) of 70 R at the position of PSR J1643-1224. Using Equation g, 
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assuming a temperature of 7000 K and using the value of E B _ V from Schlegel et al. (1998) 
to derive the dust correction factor e T = 2.79, we calculated the EM due to the ionized gas 
in Sh 2-27 towards PSR J1643-1224 to be EM = 395 cm" 6 pc. We find that n = fff = 7.6 
cm -3 , which is broadly consistent with the median value 10.6 ± 2.8 cm -3 derived from 
background RMs for Sh 2-27, within the estimated uncertainties. We also calculate a filling 
factor, / = JtStt = 0.2 for Sh 2-27, which is comparable to the value / = 0.1 that we have 
assumed in our calculations for density and magnetic field strength. 



Discussion 



4.1. Comparison with previous results 



The strengths of magnetic fields in and around Hll regions have previously been 
investigated using a wide variety of observational methods, including Zeeman splitting of 



radio recombination lines (Troland & Heiles 1977 Roberts et al. 1991 Roshi 2007) and 



Hi, OH and H 2 maser lines (Bloemhof et al. 1992; Crutcher et al. 1996 Sarma et al 



2001 Brogan & Troland 2001). Most of these spectral lines originate in relatively dense 



regions in neutral and molecular gas, therefore none of these observational probes is suited 
to examining the B vs. n relation in warm ionized gas, where the results of existing 
observations are very poorly constrained. 

Depolarization of radio synchrotron emission can be used to probe magnetic fields in 



H II regions, provided the distance to the H II region is known. Gray et al. ( 1999 ) studied 



the Faraday rotation of Galactic synchrotron emission by the H II region W4, placing an 
upper limit on the average magnetic field strength of (Bn) < 20 /^G, where angle brackets 



indicate a line-of-sight average. A similar approach was taken by Gaensler et al. (2001) 



who modeled the depolarization of Galactic radio synchrotron emission by the Hll regions 
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RCW 94 and "void 1" (which we identify as RCW 105) with n e ~ 20 cm -3 , to obtain an 
estimate of the random magnetic field strengths B r ~5 fiG (void 1) and B r = 1.2 fiG (RCW 
94). Using the fact that depolarization is negligible towards the edges of the voids were the 
gradient in RM is most severe, and assuming that the random and uniform components of 
the magnetic field are approximately equal, they determined that Bn ~ 5 fiG in void 1. For 
RCW 94, the polarization profile was modelled and the assumption that the uniform and 
random components of the field are equal was adopted, leading to the estimated value of 

Bn ~ 1.2 fiG. 



Heiles & Chu (1980) and Heiles et al. (1981) were the first to use Faraday rotation of 



background polarized point-sources to probe the line-of-sight strength and orientation of 
the magnetic fields in H II regions. They estimated the parallel component of the magnetic 



field strengths in 4 Hll regions, finding field strengths of the order Bn ~1— 10 fiG. Heiles 



et al. (1981) constrained the magnetic field strength of Sh 2-264 to the range Bn = to 



+2.6 fiG. This agrees with the magnitude and direction for the magnetic field of +2.2 ± 
1.6 fiG derived by us, within the errors, although they used an average value of n e for 
the whole Hll region that was much lower than the median value of hq that we derived. 



More recently, Mitra et al. (2003) used the RM and DM of two pulsars to estimate the 



line-of-sight component of the uniform magnetic field in the Hll region Sh 2-205, finding 
that Bn ~+5.7 fiG. In each of these studies, only one or two polarized background sources 
was available to probe each region and the foreground correction was poorly constrained. 
In our study, we have used the Faraday rotation of a large sample of polarized extragalactic 
point sources to isolate Bn in 5 Hll regions. As we have a far larger set of measurements 



than Heiles & Chu (1980); Heiles et al. (1981) or Mitra et al. (2003), the value and scatter 



of RM in our H II regions are much better characterised. As well as employing a large 
number of background sources in our study, we have also used detailed maps of Ha surface 
brightness. We are therefore able to take many independent measurements of n e and Bn 
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for each H II region. In addition, the large number of extragalactic point sources in each 
field enable us to carry out a correction for the RM contribution of the background and 
foreground to the Hll region, thus isolating the magnetic field within a discrete volume of 
interstellar space. 



A similar technique to ours was employed by Whiting et al. (2009), who described 
a 'Faraday rotation anomaly' caused by an ionised shell surrounding the Cygnus OBI 
association. Within a 2°-5° region surrounding the massive star cluster, they confirmed 
that the RMs of background radio sources change by several hundred rad m -2 . The authors 
presented a model of an ionized plasma shell driven by stellar winds, in which the magnetic 
field of the interstellar medium is enhanced by compression. The RM distribution of the 
9 polarized sources measured by Whiting et al. (2009) could not confirm or refute the 
possibility of an enhanced magnetic field in the shell, however the intensity distribution 
of Ha emission from the WHAM survey was consistent with a limb-brightened shell 



surrounding Cygnus OBI. A similar scenario was proposed by Stil et al. (2009), who 
produced 3-dimensional MHD models of a magnetized superbubble and predicted the 
distribution of RMs and the gas density that would result from this bubble. They found 
that, if the line of sight is along the Galactic magnetic field, the largest rotation measures 
are associated with a thick shell at the rim of the H II region. Comparing this model to our 



data, the RM data from Taylor et al. (2009) show no evidence for the existence of a shell 
with an enhanced magnetic field in any of our 5 H II regions, although we cannot rule out a 
thin shell. To measure an enhanced magnetic field in a thin shell would require a far denser 
grid of background RMs. We searched for evidence Ha emission from an ionized shell in 
our H II regions, but in each case we found no systematic increase in Ha surface brightness 
towards the rim. The only exception was Sh 2-220, where the Hll region is interacting 
with a molecular cloud, producing a sharp edge to the Ha surface brightness. For these 
reasons, we do not consider a magnetized shell model to reflect the data in our sample of 
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large-diameter H II regions. 



4.2. The reliability of RMs from the Taylor et al. (2009) catalogue 

Mao et al. (2010) claimed to find discrepancies between the RMs from Taylor et al. 
(2009) produced by angle fitting between two adjacent spectral channels and their own 
RM measurements of the same sources, calculated using multi-channel RM synthesis. This 
difference was attributed by Mao et al. (2010) to incorrect RM assignments in the Taylor 
catalogue caused by multiple RM components that produce non-linearity in the polarization 
angle against A 2 relation. Stil et al. (2011) addressed to these claims in Section 4.2 of their 
paper, reporting that the variance in the correlation of RMs common to the Taylor and Mao 
catalogues was underestimated compared to the errors quoted in each catalogue. Stil et al. 
(2011) also pointed out that in the Mao et al. (2010) RM sample, RMs with amplitude 
more than 20 rad m -2 are much more likely to have a large error, which introduces a 
correlation between RM error and RM amplitude. This appears to explain the discrepancy 
between the two data sets, and also indicates that the uncertainties in RMs that satisfy 
|RM| < 8 rad m~ 2 are approximately 22% larger than quoted in the Taylor catalogue. If 
this interpretation is correct, it would have the effect of increasing the smallest error bars on 
our plots by 22% (Figure 2). Given the small size of the error bars for these sources, such 



an increase would have a negligible effect on our analysis. An observational study by Law 



et al. (2011) measured the RMs of 37 polarized sources using RM synthesis on data from 
the Allen Telescope Array. They compared their results with published RMs from Taylor 
et al. (2009) and found that in that catalogue, the RMs of sources with fluxes above 200 



mJy do not suffer from nir ambiguities and have at least a 95% reliability. Van Eck et al. 



(2011) also compared their RMs to those of Taylor et al. (2009), finding a 96% correlation 



between the data sets. 
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4.3. Magnetic vs. thermal pressure in Hll regions 



The magnetic field of the Galaxy is largely confined to the disk by the pressure of the 
gas threaded by the field (Parker 1966). Measuring the individual components of the total 
pressure in the ISM is important in testing multi-phase models of hydrostatic support in 
the Galactic disk, allowing simulations of the regulation of galactic star formation rates 



(Ostriker et al. 2010) and determination of the scale height of galactic disks (Koyama & 



Ostriker 2009). Using our measurements of magnetic field strength and electron density, we 



were able to estimate the magnetic and thermal pressure in five H II regions that lie in the 
local Galactic disk. 

Figure [3] shows the relationship B\\ vs. n, for a total of 93 measurements of polarized 
sources behind 5 Hll regions. For each data point we calculated the thermal pressure, 
magnetic pressure and hence their ratio (3 th . The median P ma g, Pther, 0th an< ^ the ran g es 
that encompass 68% of the values of /3 th , centered on the median for each Hll region 



are listed in Table |3j Note that quoted j3 th = (Pther /Pmag) for all the data points in 
an Hll region, not the ratio Pther/ Pmag for that region. The black line on Figure [4] 
indicates equipartition between thermal and magnetic pressure, 2n kT e = -^ Et , assuming a 
temperature of 7000 K and a 100% ionization fraction. 

Our data reveal a wide range of values of 0th i n our 5 Hll regions, with median values 
2.0 < ]3 th < 22.4. In Sh 2-264 and Sh 2-171 ]3 th » 1 within the estimated uncertainties 
and in the other three HII regions, (3 th ~2. Values of /3th > 1 suggest that the Hll regions 
are still in a thermally-evolving/expanding phase, i.e. their morphology is dominated by 
thermal pressure and is not yet shaped by the magnetic field of the surrounding medium 



(Krumholz et al. 2007). In molecular clouds and the cold neutral medium, typical values of 



f3 t h are ~ 0.04 and ~0.3, respectively (Crutcher 1999 Heiles & Crutcher 2005). Improving 



the method by which measurements of electron density are made (e.g. by studying multiple 
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radio recombination lines or measuring radio continuum emission) may help to reduce 
the uncertainties in (5 t h- Another way to achieve this is to improve the number density 
of background polarized radio sources, thus allowing a more sophisticated removal of the 
foreground RM and reducing the measurement error in Bn. A density of 100 polarized 
sources per square degree will likely be achieved by a planned survey of the RM sky using 



the Australian Square Kilometer Array Pathfinder (Gaensler et al. 2010). 



The strengths of the magnetic fields in our H II regions are similar to estimates of the 
magnetic field strength in the diffuse ISM. This suggests that the magnetic field that we 
measure in these regions is simply the magnetic field of the warm ionized medium of the 
Galaxy. The influence of the OB star at the center of the Hll region is simply to ionize 
the surrounding few tens of parsecs, allowing us to use the Faraday rotation of background 
sources to measure the properties of the Galactic magnetic field along the line of sight. In 
this way, we have measured the strength and orientation of the Galaxy's magnetic field in 
five discrete, isolated regions of the warm ionized medium. Future work on a larger sample 
of H II regions will assist in constraining (3 t h m a wider range of H II regions. 



4.4. The magnetic field — density relation 

In Section 1 we introduced the B vs. n relation and its relevance to molecular cloud 
formation and the large-scale structure of the ISM. We now measure the slope of B vs. n 
in the low-density ISM by adding our data to existing observations of magnetic fields in 
atomic and molecular gas. 

The relationship between magnetic field and density in neutral gas and molecular 



clouds is well constrained by a wealth of observations (see, e.g. Crutcher et al. 2003). This 



relationship in the ISM is well described by B oc n K , where the exponent k has a value 
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of approximately 0.4—0.6 in molecular clouds (e.g. Crutcher 1999 ) 2 . In low density gas 



traced by neutral hydrogen and the Faraday rotation and dispersion of pulsar signals, 
existing observations suggest that the strength of the magnetic field is invariant across three 



orders of magnitude in density (Troland & Heiles 1986 Heiles & Crutcher 2005 Crutcher 



2007). However, the B vs. n relationship in the diffuse ISM is very poorly constrained by 



observations, (see Section [T]) with a sparsely sampled density range and many observations 
being upper limits. Our results from low-density Hll regions (where no ~ 0.8 — 30 cm -3 ), 
in addition to existing measurements of magnetic field strengths in the ISM allow us to 
parameterize the B vs. n relationship across this range of densities. 

Figure p| shows a plot of Bn vs. n in the density range 1 — 10 7 cm -3 , using a 



compilation of published data from Crutcher et al. (2010) (green circles and blue squares) 
overplotted with our data from extragalactic point sources seen through H II regions (red 
diamonds). The plotted density, n, is that of hydrogen atoms and molecules in neutral and 
molecular gas respectively and is the electron density n in Hll region clumps. We have 
plotted B\\, the quantity directly measured both by Zeeman splitting and Faraday rotation. 
The error bars on our data arise from the uncertainties in measurements of RM and Ha 
surface brightness. There are additional uncertainties in E#_y and /, which were shown in 



Figure pi We employed the Levenberg-Marquardt iteration technique (More 1978) to carry 



out separate parametric non-linear least squares fits to the data in the dense molecular gas 
and in the low-density (Hll + Hi) gas. Using a function of the form |B||| = An K in the 
density range n > 480 cm -3 (i.e. for the neutral and ionized medium), our fit yielded k = 
0.66 ± 0.05, which is consistent with the slope quoted by Crutcher et al.| ( |2010 ) for the same 
data. In the density range n < 480 cm" 3 , populated by Hi data and our new measurements 



2 The relationship B oc <7 Vtto tn K , where a v jot is the 1-dimensional velocity dispersion of the 



gas is also sometimes used (Basu 2000). 
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of Hll regions, our fit yielded k = 0.11 ± 0.09. Quoted uncertainties represent a oo/o 
confidence interval. 

Our data have allowed us to probe magnetic fields in discrete volumes of ionized gas 
with densities in the range 0.8 — 30 cm" 3 . The value of the density exponent k that we 
derived for the neutral plus ionized gas is very close to, but sightly above zero. This means 
that neither the atomic hydrogen in the cold neutral medium nor the ionized gas in clumps 
within our H II regions is strongly coupled to the magnetic field. It may be that there is no 
coupling (ft = 0) and we have underestimated our uncertainties. It seems more likely that 
we have chosen our density cutoff slightly higher than the limit for self-gravitating clouds 



(for example, Crutcher et al. 2010, chose a cutoff point at n= 300 cm -3 ). Our data have 
significantly expanded the density parameter space for which magnetic field strengths are 
known. Further observations of ionized gas at even lower densities will help to constrain 
this slope further. Such studies would benefit from using a measure of electron density 
that does not suffer from extinction, such as radio recombination lines or free-free radio 
continuum emission. The uncertainties in our measured values of B\\ were affected to a 
large degree by the dust extinction correction, e T , as well as the filling factor, /. Although 
the uncertainties in / are very difficult to remove, maps of radio continuum emission would 
allow a more accurate determination of the total electron density along each line of sight. 



4.5. Magnetic field orientation in the local Galaxy 



Several authors have mapped out the magnetic field of the Galaxy using measurements 



of pulsar RMs and DMs (Lyne & Smith 1989; Rand & Lyne 1994 Han et al. 2006 



Noutsos et al. 2008 Han & Zhang 2007) although the uncertainties in the distances to 



most pulsars, which depend on models of the Galactic electron density (Taylor & Cordes 



1993 Cordes & Lazio 2002), place strong constraints on the accuracy of the derived 
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magnetic field configuration. Studies of the Galactic magnetic field have also been made by 



measuring Faraday rotation of extragalactic point sources (Brown et al. 2007; Roy et al 



2008 


Taylor et al. 


2009 


Mao et al. 


2010 



electron density-weighted line-of-sight magnetic field between the pulsar and earth, whereas 
individual extragalactic RMs probe the electron density-weighted average magnetic field 
throughout the entire Galaxy. Both of these techniques are powerful and complementary, 
provided that the distances to pulsars are known well and the number of extragalactic 
point-sources is sufficient to subtract out local fluctuations in B\\ and n e . What they cannot 
do well is measure the configuration of the magnetic field in the solar neighbourhood. Here 
in H II regions, we have measured Bu in five volumes of the local Galaxy spread around three 
quadrants, thus providing a valuable addition to these existing data on the configuration of 
the local magnetic field within 1 kpc of the Sun. At the time of writing, there were also 
25 pulsars with measured RMs that have distance estimates from parallax measurements 
within the local volume. Using these combined data, we have attempted to map out the 
magnetic field configuration in the solar neighbourhood. 

Figure [5] shows a map of the local 1 kpc of the Galaxy, viewed from three orthogonal 
directions, onto which we have plotted our 5 Hll regions and all the pulsars with both 
measured RMs and with published distances measured by parallax. Both above and below 
the Galactic plane, our measurements of Bu are consistent with a preferred direction of B 



towards I ~ 90° within 1 kpc of the Sun. This agrees with the conclusions of Manchester 



(1972, 1974); Thomson fc Nelson (1980); Han fe Qiao (1994); Rand fe Lyne (1994); Han 



et al. (2006); Brown et al. (2007) and Noutsos et al. (2008), who found that the uniform 



magnetic field in the Galactic plane is oriented counterclockwise in Sagittarius-Carina spiral 
arm. Our data do not allow us to comment on the proposed reversal in the Scutum-Crux 



arm, proposed by Lyne & Smith (1989); Rand & Lyne (1994); Frick et al. (2001); Brown 



et al. (2007). Previous studies of large pulsar datasets (Han & Qiao 1994; Han et al. 1999) 
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have claimed to see a weak vertical field component in the local neighbourhood. A study 
of extragalactic point-source RMs found a preferred vertical magnetic field direction below 



the Galactic plane but no preferred direction above it (Mao et al. 2010). Our RMs (X-Z 
and Y-Z projections) show no evidence of a preferred uniform magnetic field direction 
orthogonal to the Galactic plane in the solar neighbourhood. 



5. Conclusions 

Using the RMs of an ensemble of background polarized sources and maps of surface 
brightness from published Ha surveys, we have measured the line-of-sight component of the 
magnetic field strength in 5 Galactic H II regions. Within each H II region, the line-of-sight 
component of the magnetic field is coherent. The median strength of Bu in each region 
lies between 2.2 and 6.3 /xG and the estimated median electron density ranges from 1.5 to 
14.1 cm -3 , where the electrons are assumed to reside in clumps that occupy 10% of each 
sightline through the Hll region. For one of the Hll regions, Sh 2-27, we also verified our 
estimate of B\\ and n using the RM of a background pulsar. 

With such a large number of background polarized sources, we have been able to 
subtract the RM due to the Galaxy and determine an independent measure of electron 
density and B\\ for each sightline for which we had a background RM. Our measurements 
offer a significant improvement on existing published data (or upper limits) in the 
low-density ionized ISM. Future, similar studies of Hll regions with a larger range of 
densities will constrain the relationship between Bu and no further. 

Our measurements of B\\ and uq allowed us to estimate the magnetic and thermal 
pressure in 5 Galactic Hll regions. In two regions the thermal pressure was dominant 
and in the other three, thermal and magnetic pressure were approximately equal within 
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the relatively large uncertainties of the measurements. This suggests that the H II regions 
are moving towards an advanced stage of their formative phase, where the dominance of 
thermal pressure is diminishing and the external magnetic field is able to take a role in 
shaping the regions. 

The consistent sign (positive or negative) of the RM within each H II region indicates 
that the magnetic field is largely coherent. As a consequence, we were able to measure 
the direction of B\\ in five regions of the local Galaxy. Both above and below the Galactic 
plane, the magnetic field is oriented roughly towards 1—90°. There is no preferred vertical 
direction of the Galactic magnetic field apparent from our data. Future studies involving 
measurements of B\\ in a larger number of Hn regions will extend our understanding of the 
magnetic field in the local Galaxy. 
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Table 2. Boundary conditions and details of extrinsic RM and EM subtraction for the H II regions. Pixels that lie within 
an 'On Region Boundary' (column 2) in Galactic co-ordinates that also have lua greater than 'I// a On' (defined in column 3) 
are defined as being part of the Ho region. Pixels that lie within the 'Off Region Boundary' (column 4) and have Ina less 
than 'Ihoi On' were used to subtract the Galactic background and foreground emission that does not originate from the Ho 
region. The 'Extrinsic Ina' is the mean Ho surface brightness in the off region. Where an extrinsic RM gradient is quoted, I 

is the galactic longitude, in degrees. 
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Boundary 
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40 


4° 5 < r < 6° 
I, b: +164°. 5,-12° 


17.4 


+ 18.8 


SHASSA 


Sivan 3 


r <4°.5 
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Fig. 1. — : Rotation measures (circles) and smoothed Ha surface brightness with no extinc- 
tion correction applied (grayscale) in (a) Sh 2-27, (b) Sh 2-264, (c) Sivan 3, (d) Sh 2-171, 
& (e) Sh 2-220. Maps are labelled in Galactic co-ordinates. Blue and red circles indicate 
negative and positive RMs, respectively and the diameter of the circles are proportional to 
|RM|. The legend gives 5 examples of circle sizes and the corresponding |RM| values in 
rad m~ 2 . The green circle denotes the position of PSR J1643-1224. The dashed line denotes 
the boundary of the Hll region, as defined in Table [2] The grayscale varies linearly between 
given minimum and maximum values and the colorbar shows the observed in rayleighs. 
Note that the dashed line indicating the boundary of each region is a guide, as the Ha images 
have been smoothed. 
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Fig. 1 — : continued (b) Sh 2-264 
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Fig. 1. — : continued (c) Sivan 3 
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Fig. 1 — : continued (d) Sh 2-171 
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(e) Sh 2-220 
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Fig. 1 — : continued (e) Sh 2-220 
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Fig. 2.—: l Ha vs. RM for (a) Sh 2-27, (b) Sh 2-264, (c) Sivan 3, (d) Sh 2-171 & (e) Sh 
2-220, showing the presence of coherent magnetic fields. Black points: extragalactic sources 
intersecting the sky plane outside the boundary of the H II region, red and blue circles: point 
sources that intersect the H II region and that have positive and negative RMs, respectively. 
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Fig. 2.—: continued, (b) Sh 2-264 
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Fig. 2.—: 



continued, (c) Sivan 3 
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continued, (d) Sh 2-171 
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Fig. 2.- 



continued, (e) Sh 2-220 
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Fig. 3. — : Parallel component of the magnetic field strength against electron density for sightlines 
through 5 Galactic Hn regions. Data have been corrected using the assumption that the dust 



contribution measured by Schlegel et al. (1998) lies completely in front of each Hn region. Colored 
arrows (left) indicate the displacement of any single data point if the Ha optical depth were reduced 
by a factor equal to the standard deviation of the optical depth across that H n region. Data are 
plotted assuming a filling factor of / = 0.1. The black arrow (bottom right) indicates the movement 
of each data point resulting from a change in the electron filling factor, /, by a factor of 2.5 in either 
direction. The solid black line represents equilibrium between the magnetic and thermal pressure, 
assuming T=7000 K. 
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Fig. 4. — : Bn vs. n covering nine orders of magnitude in density, using data on Hn regions, 
neutral hydrogen clouds and dense molecular cores. The Hi and molecular data are from Crutcher] 
et al. ( |2010 ) and the data from Hn regions are our data. Separate best fits to the data (black 



lines) above and below n = 480 cm 3 are shown, along with dashed lines indicating the slopes of 
la greater and less than the least squares best fit. 
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Fig. 5. — : Orientation of the magnetic field in the solar neighbourhood, shown in three orthogonal 
projections. The X-Y plane is parallel and the Z-axis is perpendicular to the Galactic plane. The 
Sun lies at the origin. Galactic longitude increases from zero (solid line) in a counterclockwise 
direction (top left panel). Blue and Red symbols indicate magnetic fields oriented away from and 
towards the earth, respectively. The square symbols represent our data for Hn regions. The 



pulsars with parallax distance measurements (Gwinn et al. 1986; Bailes et al. 1990, Kaspi et al. 



19941 |Brisken et aTl|2002| |Dodson et aL\ [20031 |Brisken et all [20031 |Ghatterjee eTe!\ [20041 [2009 



Deller et al. 2009; Verbiest et al. 2009) are represented by crosses. Circles indicate sources for 



which the rms uncertainty in Bu exceeds the magnitude of Bu. The size of each symbol represents 
the magnitude of Bu and the black lines indicate the rms uncertainty in the distances to pulsars 
measured by the annual parallax method. 



